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The Be/X-ray Transient V0332+53: Evidence for a tilt 
between the orbit and the equatorial plane? 
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ABSTRACT 

We present optical and infrared observations of BQ Cam, the optical counterpart to 
the Be/X-ray transient system V0332+53. BQ Cam is shown to be an 08-9Ve star, 
which places V0332+53 at a distance of ~ 7 kpc. Ha spectroscopy and infrared pho- 
tometry are used to discuss the evolution of the circumstellar envelope. Due to the 
low inclination of the system, parameters are strongly constrained. We find strong evi- 
dence for a tilt of the orbital plane with respect to the circumstellar disc (pressumably 
on the equatorial plane). Even though the periastron distance is only rs 10i?,*, dur- 
ing the present quiescent state the circumstellar disc does not extend to the distance 
of periastron passage. Under these conditions, X-ray emission is effectively prevented 
by centrifugal inhibition of accretion. The circumstellar disc is shown to be optically 
dense at optical and infrared wavelengths, which together with its small size, is taken 
as an indication of tidal truncation. 

Key words: binaries: general - stars: emission-line, Be - stars: individual: BQ Cam 
- pulsars:general - infrared: stars - X-rays: stars. 



1 INTRODUCTION 

The hard X-ray transient V0332+53 (X 0331+53) was first 
detected by the Vela 5B satellite during a bright outburst in 
1973 (Terrel & Priedhorsky 1984). It was rediscovered ten 
years later by Tenma during a series of smaller outbursts 
(Tanaka et al. 1983). EXOSAT observations were used to 
determine that the source pulsates with a period of 4.4 s 
(Stella et al. 1985). Doppler shifts in pulse arriving times 
indicate that the pulsar is in a 34.25-d binary orbit with an 
eccentricity e = 0.31 (Stella et al. 1985). The optical coun- 
terpart was identified by Honeycutt & Schlegel (1985) as 
the heavily-reddened early-type star BQ Cam. This object 
was observed to display highly variable Ha emission (Cor- 
bet, Charles & van der Klis 1986, and references therein) 
and infrared excess (Coe et al. 1987). These characteristics 
are typical of a Be/X-ray binary. 

In this subclass of Massive X-ray Binaries, the X-ray 
emission is believed to be due to accretion of matter from a 
Be star by a compact companion (see White, Nagase & Par- 
mar 1995; Negueruela 1998). The name "Be star" is used as 
a general term describing an early-type non-supergiant star, 
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which at some time has shown emission in the Balmer series 
lines (Slettebak 1988, for a review). Both the emission lines 
and the characteristic strong infrared excess when compared 
to normal stars of the same spectral types are attributed to 
the presence of a circumstellar disc. Most Be/X-ray binaries 
have relatively eccentric orbits and the neutron star com- 
panion is normally far away from the disc surrounding the 
Be star. Due to their different geometries and the varying 
physical conditions in the circumstellar disc, Be/X-ray bina- 
ries can present very different states of X-ray activity (Stella, 
White & Rosner 1986). In quiescence, they display persis- 
tent low- luminosity (L x 10 36 erg s _1 ) X-ray emission or 
no detectable emission at all. Occasionally, they show series 
of periodical (Type I) X-ray outbursts (L x » 10 36 - 10 37 erg 
s" 1 ), at the time of periastron passage of the neutron star 
(e.g., A 0535+26, Motch et al. 1991). More rarely, they un- 
dergo giant (Type II) X-ray outbursts (L x <; 10 37 erg s _1 ), 
which do not show clear orbital modulation. Some systems 
only display persistent emission, but most of them show out- 
bursts and are termed Be/X-ray transients. 

Like most other Be/X-ray transients, V0332+53 has 
shown both types of outbursts. The 1973 outburst lasted 
~ 100 days and peaked at ~ 1.6 Crab near July 10. It was 
clearly a Type II outburst, even though Whitlock (1989) 
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Figure 1. Composite spectrum of BQ Cam in the classification region. The short- wavelength end was taken on August 14, 1997, using 
ISIS on the WHT equipped with the R1200R grating and the Loral 1 CCD. The long-wavelength stretch was taken on November 17, 
1997, with the same instrument equipped with the R1200B grating and the EEV#10 camera. The comparison spectra (dashed) are those 
of the 09III standard i Ori (top) and the 08V standard HD 48279, both taken from the digital atlas of Walborn & Fitzpatrick (1990). 
Also displayed (bottom) is the spectrum of the 09IIe star HD 333452, from Steele et al. (1999). All spectra have been smoothed with a 
a = 0.6A Gaussian function for display. Note the strong interstellar lines in the spectrum of BQ Cam, marked with a '*'. 



found an underlying orbital modulation when the main 
trend was removed. On the other hand, the three weak out- 
bursts separated by the orbital period observed in 1983 are 
Type I outbursts. During these outbursts, the pulsed frac- 
tion was small (10 — 15%) and the temporal behaviour was 
dominated by random rapid fluctuations (Makishima et al. 
1990a). The spectrum was fitted with a power-law modified 
by cyclotron absorption. Unger et al. (1992) found that the 
pulse profile varied between a double-peaked and a single- 
peaked structure. The equivalent hydrogen column density 
remained relatively constant at ~ 1 x 10 22 atom cm -2 . A 
prominent absorption feature at 28.5 keV, if attributed to 
electron cyclotron resonance, implies a magnetic field at the 
surface of the neutron star of 2.5 x 10 12 G (Makishima et al. 
1990b). 

A new outburst was discovered by Ginga in September- 
October 1989. The source remained very bright for more 
than two weeks, indicating that this was a Type II out- 
burst. A quasi-periodic oscillation, possibly implying the 
presence of an accretion disc, was detected (Takeshima et 
al. 1994). V0332+53 has not been detected by the BATSE 
experiment since the CGRO satellite started operations in 
April 1991 (Bildsten et al. 1997). It is not detected with any 
significance by the All-Sky Monitor (ASM) on board RXTE 
either, according to the quick-look results provided by the 
RXTE I ASM team. Stella et al. (1986) interpret the lack of 
quiescence emission as an effect of centrifugal inhibition of 
accretion. 



2 OBSERVATIONS 

We present data obtained as a part of the Southamp- 
ton/Valencia/SAAO long-term monitoring campaign of 
Be/X-ray binaries (see Reig et al. 1997a), consisting of opti- 
cal spectroscopy, infrared and optical broad-band photome- 
try of BQ Cam, the optical counterpart to V0332+53. 



2.1 Blue optical spectroscopy 

The source was observed on August 14, 1997 using the 4.2- 
m William Herschel Telescope (WHT), located at the Ob- 
servatorio del Roque de los Muchachos, La Palma, Spain. 
The blue arm was equipped with the Lorall CCD and the 
R1200R grating, which gives a nominal dispersion of ~ 0.25 
A/pixel. A second observation was taken on November 14, 
1997. On this occasion the blue arm was equipped with the 
R1200B grating and the EEV#10 CCD, giving a nominal 
dispersion of ~ 0.22 A/pixel over ~ 900 A. A composite 
spectrum is shown in Fig. |l| The signal-to-noise ratio (SNR) 
of the August spectrum is relatively low. The Loral camera 
introduces several artifacts in the range A 4150-4250 A, 
where the spectrum could well be dominated by noise. A 
very strong spurious feature was present at the wavelength 
where He II A 4200 A should be found. Between A 4320 A and 
A 4500 A, where both spectra overlap, all the features look 
very similar and only the higher quality November spectrum 
is shown. 
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Table 1. Details of the Ha spectroscopy. The FWHM of Ha has 
been calculated by fitting a single Gaussian to the line profile and 
should therefore be taken only as an approximation. In the higher 
resolution spectra, however, this procedure gives values consistent 
with direct measurement to within 5%. 



Date of Telescope Dispersion EW of FWHM 

Observation(s) (A/pixel) Ha (A) (km s _1 ) 



Jan 28, 1990 


INT 


~ 0.2 


-4.5 


250 


Sep 02, 1990 


INT 


~ 0.8 


-5.0 


250 


Nov 14, 1990 


INT 


~ 0.8 


-5.7 


260 


Dec 27, 1990 


INT 


~ 0.4 


-6.8 


260 


Jan 27, 1991 


INT 


~ 0.8 


-6.5 


240 


Aug 28, 1991 


INT 


~ 0.4 


-4.6 


280 


Dec 14, 1991 


INT 


~ 0.4 


-4.8 


270 


Aug 18, 1992 


PAL 


~ 0.9 


-3.1 


230 


Nov 13, 1992 


PAL 


~ 0.9 


-3.8 


290 


Dec 05, 1993 


PAL 


~ 0.9 


-6.8 


280 


Aug 14, 1997 


WHT 


~ 0.4 


-4.2 


280 


Nov 14, 1997 


WHT 


~ 0.4 


-3.9 


270 



2.2 Red optical spectroscopy 

We have monitored the source since 1990, using the 2.5- 
m Isaac Newton Telescope (INT) and 4.2-m WHT, both 
located at the Observatorio del Roque de los Muchachos, 
and the 1.5-m telescope at Palomar Mountain (PAL). A log 
of observations, together with some parameters measured, 
is presented in Table ^ The Palomar spectra have relatively 
low SNR and the error in EW, arising due to the difficulty of 
determining the continuum, is ~ 15%. The line shapes are 
also difficult to establish. In contrast, the INT and WHT 
spectra have all relatively high resolution and errors in the 
EW of Ha are ~ 5%. The double-peaked structure of Ha 
is only clearly visible in the spectra with dispersions of 0.4 
A/pixel or better. 

All the data have been reduced using the Starlink soft- 
ware packages CCDPACK (Draper 1998) and FIGARO (Short- 
ridge et al. 1997) and analysed using FIGARO and DIPSO 
(Howarth et al. 1997). 

Ha spectra normally show He I A6678 A as a very weak 
emission feature when the SNR is sufficiently large for it 
to be separated from the noise. An analysis of Ha variabil- 
ity in BQ Cam during 1990-1991 has been presented in 
Negueruela et al. (1998). Many of the spectra whose param- 
eters are listed in Table [j] are shown in their Fig. 8, and 
are therefore not reproduced here. Negueruela et al. (1998) 
found that during 1990-1991, the Ha line presented V/R 
variability with a quasi-period of ~ 1 year, but this variabil- 
ity stopped late in 1991. 

2.3 Optical Photometry 

Optical photometry of the source was obtained on Novem- 
ber 11, 1997, using the 1-m Jakobus Kapteyn Telescope 
(JKT) at the Observatorio del Roque de los Muchachos, 
La Palma, Spain. The telescope was equipped with the Tek4 
CCD and the Harris filter set. Conditions were photometric. 
Instrumental magnitudes were extracted through synthetic 
aperture routines contained in the iraf package, and trans- 



formed to the Johnson/Cousins system through calibrations 
derived from observations of a number of Landolt (1992) 
standard stars taken on the same night. The values measured 
are U = 17.74 ± 0.10, B = 17.29 ± 0.02, V = 15.73 ± 0.02, 
R = 14.69±0.02 and I = 13.27±0.08. Errors in the U and I 
bands are dominated by calibration uncertainties (mostly in 
the colour correction equations). The smaller errors in B,V 
and R are dominated by measurement errors. 

2.4 Infrared Photometry 

Infrared observations of BQ Cam are listed in Table ^. They 
have been obtained with the Continuously Variable Filter 
(CVF) on the 1.5-m Carlos Sanchez Telescope (TCS) at the 
Teide Observatory, Tenerife, Spain, and the UKT9 detector 
at the 3.9-m UK Infrared Telescope (UKIRT) at the Mauna 
Kea observatory on Hawaii. The December 1994 observation 
was taken with IRCAM mounted on UKIRT and 540-s total 
exposure in each filter. Data for the period 1983 - 1986 have 
already been reported in Coe et al. (1987). The observations 
presented here extend for a further ten years. The long-term 
lightcurve is shown in Fig. @. 



3 RESULTS 

3.1 Spectral classification 

The spectrum of BQ Cam in the classification region is dis- 
played in Fig. |IJ We can see that H/3 is in emission, with 
a clear double-peaked structure (see also Fig. almost 
to the continuum level. The He I lines at AA 4713 & 5016 
A also show weak double-peaked emission, just above the 
continuum level. In contrast, H7 only shows two very weak 
emission components on the wings of the absorption line. 

The very strong diffuse interstellar lines are consistent 
with the high reddening of the object. The strength of the 
He 11 lines clearly identifies the object as an O-type star. This 
is confirmed by the strong Si IV lines. Overall the spectrum 
is very similar to that of LS 437, the optical counterpart to 
3A 0726-260 (Negueruela et al. 1996). 

In Fig. [j], the spectrum of BQ Cam is compared to 
that of two MK standards from the digital atlas of Wal- 
born & Fitzpatrick (1990), t Orionis (09III) and HD 48279 
(08V). Also shown is a high SNR spectrum of the Be star 
HD 333452 (09IIe) from Steele, Negueruela & Clark (1999). 
Points to be noted are: 

• The H 1 and He I lines are weaker in BQ Cam than in 
any standard, due to the presence of circumstellar emission. 

• BQ Cam must be later than 07, since Hei A4471 A>> 
Hen A4541 A. 

• The presence of He 1 AA 4026, 4144 A the strength of 
the Si iv doublet, Cm A4072 A and Hei A4388 A all point 
to BQ Cam being no earlier than 08. 

• The main luminosity criterion, namely, the strength of 
the Si IV AA 4089, 4116 A doublet is difficult to judge in 
BQ Cam, because Hei A4121 A is not present (presumably 
filled-in by emission) and the quality of the spectrum is low 
at that end. 

• The Mgn A4481 A line is visible at this resolution in 
both evolved 09 stars, but it is absent in BQ Cam, indicat- 
ing a lower luminosity class. Similarly, Nm AA 4511-4515 
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Figure 2. Infrared lightcurves of BQ Cam during 1983-1995. J magnitudes are represented by open squares, H by stars, K by filled 
squares and L' by triangles. Error bars have been removed for clarity (see Table |^ for details). 



A is much weaker in BQ Cam than in the giants. This line 
is also stronger in HD 48279, but Walborn & Fitzpatrick 
(1990) warn that this object is N-enhanced. 

• The N ill lines are only seen in absorption in the 08 - 09 
range in main-sequence stars (Walborn & Fitzpatrick 1990) . 
Nm AA 4379, 4642 A and possibly Nm A4511 A are visible 
on the spectrum of BQ Cam. 

• The equivalent width (EW) of He II A4686 A in BQ Cam 
(0.6 A) is inside the range typical for 08-09 stars (Conti 
& Alschuler 1971). 

• The O II A4367 A line, which can be seen in the spectra 
of the two 09 stars is also visible in the spectra of the 09V 
standards 10 Lac and HD 93028 at this resolution, hinting 
that BQ Cam could be earlier. 

• The complete absence of C ill A4650 A in BQ Cam is 
surprising. In main sequence stars, it is clearly visible as 
early as 08 (Walborn & Fitzpatrick 1990). This line is also 
absent in LS 437. There is no reason to think that BQ Cam 
is carbon deficient, since C ill A4072 A is present and there is 
no sign of N-enhancement. One possible explanation is that 
the inner regions of the circumstellar discs of Oe stars are 
hot enough to produce C ill emission, which would be filling 
in the A4650 A line. 

From all the above, it is clear that BQ Cam is an un- 



evolved star in the 08-9 range. An 08.5V classification 
is very likely, but, since the presence of emission affects the 
main classification criteria, we prefer to be cautious and give 
an spectral type 08-9Ve. 



3.2 Distance 

In order to obtain an estimate of the distance to the system, 
it is necessary to determine the interstellar absorption in 
its direction. This is complicated because Be stars present 
circumstellar reddening, due to the infrared continuum emis- 
sion. The calculation by Kodaira et al. (1985) of Ay = 7.4 
is a gross overestimate, since they took the infrared mag- 
nitudes measured in 1983 to be those intrinsic to the star, 
while our data show that the disc emission accounted for at 
least ~ 0.8 mag in J. Our UBVRI photometry was taken 
at a time when the circumstellar disc was certainly small 
(see Sections 3.3 and|^). The small EW of Ha at the time of 
the observations (see Table should be accompanied by a 
small circumstellar reddening (see Dachs, Engels & Kiehling 
1988). 

Since the intrinsic colour of a late O-type star is 
(B — V)o = —0.31 for luminosity classes III-V (Schmidt- 
Kaler 1982), the measured (B — V) = 1.56 ± 0.03 implies 
E(B — V) — 1.87. This value is almost identical to the 
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Date 




TJD 


J 




H 






K 




L' 


1983 


Nov 23 


4f)662fU1 


1 1 S9+0 03 


10 


70+0 


03 


10 


31+0 


03 


9 6+0 05 


1983 


Dec 23 


45RQ2CU1 


1 1 S9+0 03 


10 


60+0 


03 


10 


20+0 


03 




1984 


Jul 22 


45904 ("TH 


1 1 44+n 03 

X X i 11 1 KJ . VJ > J 








10 


34+0 


03 




1985 


Jul 22 


46269(jn 

Tt VJ £4 VJ %J \ \J J 


1 1 fin+n 03 

X -L . V/ VJ 1_ \J . VJ ' J 


10 


94+0 


03 


10 


53+0 


03 




1986 


Oct 10 


46714(U) 


1 1 634-0 03 

X X . VJ • J _l_ \J . LJO 


11 


01+0 


03 


10 


60+0 


03 


1 n 1 q+o 05 

1U.1 <J _1_ y-J ■ \ ><J 


1986 


Oct 23 


46727fin 

4:U i i 1 I VJ ) 


1 1 65+0 fl3 

X X . \J ij _l_ KJ . L/iJ 


11 


00+0 


03 


10 


62+0 


03 


1 n 22+0 05 

XL/. KJ . \J\J 


1987 


Feb 25 


46852fin 


1 1 66+0 03 


11 


06+0 


03 


10 


64+0 


03 


1 n 21 +0 05 


1987 


Aug 18 


47026('U1 

1 VJ V/ \ VJ J 


1 1 fin+n 03 


10 


99+0 


03 


10 


62+0 


03 


9 97+0 09 


1987 


Nov 05 

1 N VJ V VJ \J 


471 05('TI1 


1 1 R4+0 03 


11 


02+0 


03 


10 


60+0 


03 


1 n 1 1 +0 1 1 


1988 


Jan 03 


47164(T) 


1 1 62+0 04 


11 


04+0 


04 


10 


74+0 


04 




1988 


Jul 23 


47366fin 


1 1 44+0 02 


10 


78+0 


02 


10 


42+0 


03 


1 n 02+0 08 

XL/. L/ii_l_L/. L/O 


1989 


Dec 15 


47876("U"l 


1 1 47+0 01 


10 


84+0 


01 


10 


46+0 


01 


10 09+0 08 


1989 


Dec 21 


47882(Tn 


11.42+0.01 


10 


80+0 


01 


10 


42+0 


01 




1991 


Jan 22 


48279(T) 


1 1 Q5+0 06 


11 


30+0 


03 


10 


84+0 


03 




1991 


Aug 23 


48492(T) 


1 1 074-n 04 

X X . <y 1 _l_L/.Wrt 


11 


40+0 


03 


10 


97+0 


03 




1991 


Aug 25 


48494(T) 


12 08+0 06 


11 


36+0 


03 


11 


01+0 


01 




1991 


Aug 27 


48496f'T > ) 


1 1 sn+0 05 

X X . OU_l_U. v/' J 


11 


28+0 


04 


10 


96+0 


04 




1991 


Nov 04. 

1 \ VJ V L/rt 


48565fin 


1 2 00+0 02 

x^i . wo ii . uzi 


11 


37+0 


01 


10 


95+0 


02 


1 n SQ+0 06 

XL/. *J Z> -J— V/ . VJ VJ 


1991 


Dec 01 


48592(T) 


12.2 ±0.2 


11 


47+0 


07 


11 


08+0 


06 




1992 


Jan 14 


48636(U) 


12.05+0.01 


11 


46+0 


01 


11 


06+0 


01 


10.44+0.04 


1992 


Feb 12 


48665(T) 


12.16+0.04 


11 


67+0 


05 


11 


40+0 


06 




1992 


Aug 17 


48852(U) 


11.99+0.01 


11 


42+0 


02 


11 


00+0 


01 


10.56+0.09 


1993 


Jan 15 


49003(T) 


11.92+0.07 


11 


15+0 


04 


10 


73+0 


03 




1994 


Nov 07 


49664(T) 


12.16+0.06 


11 


61+0 


03 


11 


21+0 


05 




1994 


Dec 14 


49701 (U) 


11.89+0.02 


11 


43+0 


02 


11 


12+0 


02 




1995 


Jan 02 


49720(T) 


11.97+0.14 


11 


28+0 


05 


11 


04+0 


06 





Table 2. Observational details of the IR, photometry. Observations marked T are from the TCS, while those marked U are from UKIRT. 
The first five observations are taken from Coe et al. (1987). 

E(B-V) = 1.88 + 0.1 deduced from the strength of different 
interstellar bands by Corbet et al. (1986), though the photo- 
metric determination is more reliable. It is also compatible 
within the errors with the value of JVh derived by Unger et 
al. (1992) from EXOSAT X-ray data taken during the 1983 
Type I outbursts, which implies E(B — V) « 1.7 (Bohlin, 
Savage & Drake 1978). Since a very small amount of cir- 
cumstellar reddening could be present, the derived distance 
should be taken as a lower limit. Zorec & Briot (1991) and 
Fabregat & Torrejon (1998) have noted that Be stars are on 
average 0.3 mag brighter than main-sequence objects (due 
to the added luminosity of the circumstellar disc). Though 
the disc surrounding BQ Cam must be small, some contribu- 
tion to the absolute luminosity should be expected. However, 
we will use the absolute magnitude of a normal 08.5V star 
My = —4.5 (Vacca, Garmany & Shull 1996), once again 
taking into account that the distance calculated will be a 
lower limit. 

Using the standard reddening of R — 3.1, we derive d — 
7.6 kpc. However, using Schmidt-Kaler's (1982) expression 
for the reddening to O stars, we find R — 3.3, which gives 
a distance of d — 6.3 kpc. Given the uncertainty in R and 
taking into account the above considerations, we will accept 
the distance d = 6 kpc as a lower limit (unless the reddening 
in that direction is exceptionally strong). An estimate of the 
different factors mentioned above, would indicate a range 
6 < d < 9 kpc for the distance to BQ Cam. 

We note that for an 09III star, Mv = —5.5 (Vacca 
© 0000 RAS, MNRAS 000, 000-000 



et al. 1996) and the implied distance is d ^ 10 kpc, which 
would place the object well outside the galactic disc. This 
is taken as confirmation of the main-sequence classification 
for BQ Cam. 



3.3 System parameters 

The very low mass function of V0332+53, f(M) = 0.10 ± 
0.03 (Stella et al. 1985) indicates that the orbit of the neu- 
tron star is seen under a very small inclination angle. As- 
suming a lower limit for the mass of an O-type compan- 
ion M» <; 20Mq and the standard mass for a neutron star 
M x = 1.44 Mq, an inclination angle i <, 10° .3 ± 0°.9 is ob- 
tained. Waters et al. (1989) have argued that the orbital 
plane of Be/X-ray binaries with close orbits should not be 
very inclined with respect to the equatorial plane in which 
the circumstellar disc is supposed to form. 

We have used our high-resolution spectra from Novem- 
ber 14, 1997 to measure the parameters of several emission 
and absorption lines, which are listed in Table ^, in order 
to estimate the rotational velocity of BQ Cam. This is par- 
ticularly difficult given the presence of a circumstellar com- 
ponent. We have selected the two strongest He lines in the 
blue. He II A4686 A is not likely to be contaminated by any 
emission component from the disc, but can be affected by 
non-LTE effects. On the other hand, He I A4471 A is likely 
to be affected by circumstellar emission, which will reduce 
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Figure 3. Emission lines in BQ Cam observed on November 14, 
1997 with the WHT equipped with the the R1200B grating and 
the EEV#10 CCD on the blue arm and the R1200R grating and 
the Tek5 CCD on the red arm. From bottom to top, H/3, He I 
A6678A (dotted line) and Ha. All the lines have been normalised 
to the continuum level (note that H/3 emission is inside the wings 
of the photospheric absorption and therefore below the continuum 
level). For the line parameters, see Table [j| Also plotted is the 
Ha profile from December 12, 1991, displaced by 0.1 units of the 
continuum. Note the similarity in shape, intensity and FWHM. 



its FWHM. Therefore any estimation of vsini based on this 
line should be taken as a lower limit. 

From Buscombe's (1969) approximation, the mea- 
sured FWHMs corrected for instrumental broadening im- 
ply vsini = 160 km s -1 from the Hell line and vsini = 
130 km s~ 1 from the He I line. Similarly, using the correlation 
between FWHM of He I A4471 A and v sin t from Slettebak 
et al. (1975), we obtain vsini — 135 km s" 1 . The two values 
derived from the He I A4471 A line are very similar and set 
a lower limit for vsini. 

A very different way of estimating the rotational veloc- 
ity is by using the mean relation between FWHM of the Ha 
emission line and vsini for Be stars from Hanuschik, Kozok 
& Kaiser (1988). Since the FWHM of Ha in BQ Cam had 
not changed significantly during 6 years, we deduce that the 
disc is dynamically stable and the correlation can be trusted 
to a relatively high degree. We note that the scatter in the 
correlation is due to the inclusion of measurements for stars 
with dynamically unstable discs. 

Using 



lot 



FWHM 



-0.2 log W a + 0.11 



2(v sin i) 

we obtain vsini = 140kms _1 . Similarly, using the mean 
relation between peak separation of Ha and vsini for Be 
stars (Hanuschik et al. 1988), 

AVpeak 



Io. 



-0.4 log W a - 0.1 



2(v sin i) 

we obtain vsini = 170 km s -1 

All the above estimates provide similar values. The esti- 



mates based on He I A4471 A yield a lower limit for v sin i 
130 km s . Averaging the four estimates gives a value of 
~ 150 km s -1 . The errors associated with this estimation 
are formally large. However, the shape and FWHM of Ha 
are not compatible with a very low vsini (^ lOOkms -1 ), 
while a value approaching v sin i ~ 200 km s _1 does not seem 
compatible with the comparatively small width (both at the 
base and at half-maximum) of Ha when a large population 
of Be stars are considered (see, for example, Hanuschik et 
al. 1996). 

Since all Be stars are believed to be fast rotators, this 
value of vsini confirms that the star is seen under a small 
inclination angle. We note, however, that in order to show 
a vsini ~ 150kms _1 with an inclination angle i = 10°, 
the rotational velocity of the star should be v ~ 900 km 
s" 1 , well above the break-up velocity of a late O-type star 
(<; 600 km s _1 ). Assuming an upper limit for the rotational 
velocity of v = 0.8vb rca k & 480 km s -1 , still gives i <; 19°. 
An orbital inclination of i = 19° , would imply an enormously 
undermassive primary with M* = 5A/q (the errors in the 
orbital parameters allow up to M* = 7Mq within 2-a). Even 
if we assume that the O star is rotating at break-up velocity, 
i = 15° implies M* = 8A/q, still undermassive by a factor 
> 2. We note the uncertainty in our estimate of vsini, but 
the constraints v < 600 km s _1 and i 10° can only be 
compatible with vsini 100 km s _1 . This is not only far 
away from our estimate of vsini, but also very difficult to 
reconcile with the FWHM and shape of Ha (see Hummel 
1994, Hanuschik et al. 1996). 

On the other hand, we have no strong reasons to ex- 
pect a very undermassive optical star. The calculations by 
Vanbeveren & De Loore (1994) show that, under certain 
circumstances, mass transfer in massive binaries can lead 
to the formation of overluminous post-main-sequence stars 
with compact companions (e.g., Vela X-l). It is not clear how 
noticeable this effect will be as long as the star which has 
received mass remains in the main sequence, and whether 
this will affect T c g (and therefore, spectral class). Gies et 
al. (1998) have found evidence suggesting that the optical 
component of the Be + sdO binary <f> Per is moderately un- 
dermassive, but no evidence exists for any main-sequence 
component of an X-ray binary being undermassive by a fac- 
tor > 2. As a consequence, we believe that the discrepancy 
in the two values for sin i strongly suggests that the orbital 
plane is not exactly aligned with the equatorial plane of the 
Be star, even though the difference may be small (~ 10°). 

Hummel (1994) has shown that, for inclination an- 
gles i 30°, the profile of emission lines from Be stars 
is dominated by the flank inflections generated by non- 
coherent scattering, giving rise to what is known as the 
wine-bottle shape. Wine-bottle shapes are found for Be stars 
with vsini 250 km s _1 and Hanuschik et al. (1996) es- 
timate that flank inflections are visible for inclinations up 
to i ~ 60° . However, the 1997 Ha profiles of BQ Cam have 
no sign of flank inflections (see Fig. ^|). Since it is not rea- 
sonable to suppose that i > 60° for BQ Cam, we interpret 
the absence of flank inflections as proof that the envelope of 
BQ Cam is small and the optical depth in the vertical direc- 
tion is not large enough to produce the wine-bottle profile 
typical of non-coherent scattering. Therefore the peak sep- 
aration of emission lines in November 1997 will reflect the 
actual extent of the disc. 
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Line 



EW (A) FWHM 
(km s" 1 ) 



Peak separation 
(km s" 1 ) 



Ha 

Hj8 

He I A6678A 
Hen A4686A 
He I A4471A 



-3.9 
- 0.8 
-0.4 
0.6 
0.6 



270 



270 
220 



155 
200 
230 



Table 3. Details of several lines in the spectrum of BQ Cam, 
observed on November 14, 1997, with the WHT. Emission lines 
are shown in Fig. M. The EW of H/3 is only approximate, due to 
the contamination by the interstellar band at A 4885 A. Note that 
He II A 4686A and He I A447lA are photospheric absorption lines, 
though the latter could be affected by a circumstellar emission 
component. 



3.4 Disc evolution 

lye & Kodaira (1985) and Corbet et al. (1986) describe 
radial-velocity changes in the Ha emission line and inves- 
tigate their possible connection with the orbital period. Our 
Ha spectroscopy shows that these changes were also present 
during 1990 - 1991, but Negueruela et al. (1998) have shown 
that these velocity variations can be explained by quasi- 
cyclic V/R variability with a quasi-period ~ 1 year. Similar 
cyclic variability is seen in many other Be/X-ray binaries 
(Negueruela et al. 1998). It is noteworthy that the system 
was displaying V/R variability in 1983 - 1984 when it showed 
a short span of X-ray activity and again in 1990, immedi- 
ately after the 1989 Type II outbursts. The possibility of a 
causal connection between V/R variability and X-ray activ- 
ity has been discussed in Negueruela et al. (1998). 

The infrared lightcurves (see Fig. ^) of BQ Cam show a 
general fading trend, only interrupted by a brief brightening 
in 1988- 1989. The brightest magnitudes observed are those 
from late 1983, when the source was active in the X-rays, 
coinciding with the highest EWs of Ha reported (—8 and 
-10 A; Kodaira et al. 1985; Stocke et al. 1985). The decline 
of the strength of Ha (lye & Kodaira 1985) was accompanied 
by the fading of infrared magnitudes. Corbet et al. (1986) 
and Coe et al. (1987) interpret the decline as being due to 
the dispersion of the circumstellar disc of the Be star. 

The infrared magnitudes remained stable during 1985 - 
1987, but brightened again in 1988, reaching values similar 
to those of 1983. After the type II outburst in late 1989, 
the infrared magnitudes faded to a deeper minimum, where 
they have remained until 1995, though showing considerable 
short-time variability. There does not seem to be any corre- 
sponding systematic change in Ha EW (see Table W). The 
variability in Ha during 1990-1991 can be associated with 
the V/R cycle seen at the time (Negueruela et al. 1998). 
Given the similarity between the relatively high-resolution 
spectra of 28 August 1991, 14 December 1991, 14 August 
1997 and 14 November 1997, it seems unlikely that any sig- 
nificant V/R variability has been present after 1991. 

It is noteworthy that, while the infrared colours have 
experimented large fluctuations (~ 0.8 mag in J and ~ 1.2 in 
K), the associated colours have remained much more stable 
(with (J — K) ~ 0.8 — 1.2). There is no clear correlation 
between the brightness and the colours. If, for instance, we 



compare the J magnitude with (J — K), we see that the 
faintest observations can be either very blue (TJD 48665) 
or relatively red (TJD 48494). The brightest observations 
are on average relatively red, but not more than some faint 
points. 

Using the correlations of Rieke & Lebofsky (1985), from 
the observed Av = R X E(B — V) = 6.17 we deduce 
Aj — 1.74 and Ak = 0.69, implying an interstellar redden- 
ing E(J — K) — 1.05. Using the fainter and bluer infrared 
observations (TJD 48665), we find J = J - Aj = 10.42 and 
K = K - A K = 10.71, implying (J - K) = -0.29 ± 0.07, 
which is roughly compatible with Wegner's (1994) average 
value of (J — K)o = —0.18 for 09V stars, if we consider the 
errors in his value and in the standard relations used. The 
measured E(J - K) = (J - K) - (J - K)o = 0.94 is, within 
the errors, compatible with the value for interstellar redden- 
ing found above and implies that no circumstellar reddening 
was present. This corresponds to a state in which the disc 
is optically thin at all infrared wavelengths and very little 
infrared emission is produced (see Dougherty et al. 1994). 
Brighter magnitudes with a similarly blue colour (as in TJD 
49701) must represent a state in which the disc is producing 
a significant amount of infrared emission, but still remains 
optically thin, giving rise to no circumstellar reddening - a 
condition which could be associated with a very small disc 
(Dougherty et al. 1994). Very faint magnitudes with redder 
colours (as in TJD 49664) represent states in which little 
emission is present, but the disc has become (partially) op- 
tically thick in K, which can be associated with a larger disc 
or with a change in the density gradient. When the disc is 
very bright (as in 1983 or 1988), it is always relatively red. 
The circumstellar emission is very intense (contributing ~ 1 
mag), but the disc is optically thick in all wavelengths, and 
the circumstellar reddening remains constant at a low value, 
presumably because the disc is still small (Dougherty et al. 
1994). 



4 DISCUSSION 

We have shown that the optical counterpart to V0332+53 
is an 08-9Ve star, further skewing the spectral distribution 
of Be/X-ray binary transients towards earlier spectral types 
(see Negueruela 1998) . Our distance estimate is higher than 
those of previous authors, who observed brighter magnitudes 
due to larger disc contamination. Honeycutt & Schlegel 
(1985) report B = 17.04 ±0.06 and (B — V) = 1.62 ±0.06 in 
two separate observations taken on November 27-28, 1983, 
and February 21, 1984, while our observations show magni- 
tudes fainter by ~ 0.3 mag, though the values for (B — V) 
are compatible within the observational errors. Normally ob- 
served photometric variability of classical Be stars is typi- 
cally ^ 0.2 mag. Therefore BQ Cam does not show vari- 
ations as large as V635 Cas, the optical counterpart to 
4U0115±634 (Negueruela et al. 1997), but it is more vari- 
able than most classical Be stars. At a distance of ~ 7 kpc, 
the maximum X-ray luminosity of V0332±53 observed by 
Vela 5B is L x J> 10 38 erg s _1 , close to the Eddington lumi- 
nosity for a neutron star. At this distance, BQ Cam is not 
likely to be part of the Perseus arm, but rather of an outer 
galactic arm (see Kimeswenger & Weinberger 1989). 

With an inclination angle i 10°, the orbital solution 



© 0000 RAS, MNRAS 000, 000-000 



8 /. Negueruela et al. 



for V0332+53 implies a x £ 8.5 x 10 10 m. For a companion 
mass M* i; 20Mq , this implies a periastron distance a pcl <; 
6.3 x 10 10 m ~ 107?*, where the radius of the optical star 
is assumed to be 7?* = 97?q (Vacca et al. 1996). Using our 
value for wsini and Huang's (1972) law 



7?d 



( 2v sin i \ 



where u pca k is the separation between the peaks of an emis- 
sion line and Keplerian rotation of the envelope is assumed 
(which gives upper limits), we obtain outer emission radii 
R d = 4.0,2.5 and 1.87?* for Ha, H/3 and He I A6678A in 
November 1997. This clearly shows that the neutron star 
does not come close to the dense regions of the circumstel- 
lar disc in the present situation. The neutron star is only 
reached by the low-density outer envelope and accretion is 
centrifugally inhibited (Stella et al. 1986). The reduced size 
of the circumstellar disc strongly points at the possibility 
of disc truncation by the neutron star, an idea advanced by 
Okazaki (1998) and supported by the results of Reig, Fab- 
regat & Coe (1997b). 

This truncation would explain why we observe instances 
of a small disc which is optically thick at all wavelengths, 
implying a very high density. We note that V635 Cas shows 
large brightness variations with little change in the associ- 
ated colours (Negueruela et al. 1997), a behaviour that could 
be associated with a very optically thick disc (Dougherty et 
al. 1994) . This variability extends to the B band, which can 
change by at least 0.6 mag. It seems likely then that the disc 
in BQ Cam is not so optically thick as that in V635 Cas at 
optical wavelengths, but can be very optically thick in the 
infrared. Roche et al. (1999) find that optical and infrared 
observations of the Be/X-ray transient Cep X-4 are also best 
explained with a truncated dense disc. 

We have established that the orbit of the neutron star 
is likely to be inclined with respect to the equatorial plane 
of the Be star (in which the disc is suppose to lie while it is 
dynamically stable). This is in contradiction with the gen- 
eral argument presented by Waters et al. (1989), but it is in 
not unexpected. We note that the Be + neutron star sys- 
tem PSR B1259— 63, which is believed to be a representative 
of the class of systems which will evolve into Be/X-ray bi- 
naries, is likely to have a tilted orbit with respect to the 
equatorial plane of the Be star (see Ball et al. 1999 and 
references therein), and the B + neutron star system PSR 
J0045— 7319, which must have formed in a way analogous to 
Be/X-ray binaries, has been shown to have a rotation axis 
misaligned with the orbit (Kaspi et al. 1996). The relevance 
of this misalignment to the formation and evolution of bi- 
nary systems containing neutron stars has been discussed by 
van den Heuvel & van Paradijs (1997) and Iben & Tutukov 
(1998) 

Corbet & Peele (1997) have suggested a possible 34.5-d 
period for the Be/X-ray binary 3A 0726— 260. If this period 
was to be confirmed, the comparison between V0332+53 
and 3 A 0726—260 would be most interesting, since the two 
systems would have neutron stars orbiting Oe stars of al- 
most identical spectral types with extremely similar orbital 
periods. In contrast, the X-ray activity of these two sys- 
tems is completely different. V0332+53 is a transient, which 
spends most time in quiescence and shows very bright out- 
bursts, while 3A 0726— 260 seems to be a persistent low- 



luminosity source with small outbursts. The difference in 
pulse periods (4.4-s against 103. 2-s) could be reflecting the 
very different behaviour of accreted material in magnetic 
fields of very different intensity. However, we note that the 
quiescent luminosity of 3A 0726—260 is almost identical to 
that of A 0535+262, which has almost the same spin period, 
but a much broader orbit. This, together with the fact that 
the source lies nowhere close to the P rb/Ps relationship for 
Be/X-ray binaries, casts some doubt on the orbital period 
until it can be confirmed using orbital Doppler shift in the 
arrival time of pulses. 



5 CONCLUSIONS 

We have presented long-term photometry and spectroscopy 
of the optical component of the Be/X-ray binary V0332+53, 
which indicate that it is an 08-9Ve star at a distance of ~ 7 
kpc. We find evidence for a tilt of the orbital plane with 
respect to the equatorial plane. The lack of recent X-ray 
activity is explained by the fact that the dense regions of 
the circumstellar disc around the Oe star do not reach the 
orbit of the neutron star. The low inclination of the orbit 
allows us to determine a periastron distance a pcr ^ 107?*, 
while measurements from our high-resolution spectroscopy 
of emission lines set the outer radius of the Ha emitting 
region at 7?d & 47?*. Under these conditions, centrifugal in- 
hibition of accretion effectively prevents any X-ray emission. 
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